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Abstract 

We analyze a simple opinion formation model consisting of two parties, A and B, and a group /, of 
undecided agents. We assume that the supporters of parties A and B do not interact among them, 
but only interact through the group /, and that there is a nonzero probability of a spontaneous 
change of opinion (A £5 /, B £5 I). From the master equation, and via van Kampen's f2-expansion 
approach, we have obtained the "macroscopic" evolution equation, as well as the Fokker-Planck 
equation governing the fluctuations around the deterministic behavior. Within the same approach, 
we have also obtained information about the typical relaxation behavior of small perturbations. 
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I. INTRODUCTION 



The last few years have witnessed a growing interest among theoretical physicists in 
complex phenomena in fields departing from the classical mainstream of physics research. 
In particular, the application of statistical physics methods to social phenomena has been 
discussed in several reviews jl|, 0, S, U, 0] . Among these sociological problems, one that has 
attracted much attention was the building (or the lack) of consensus. There are many differ- 
ent models that simulate and analyze the dynamics of such processes in opinion formation, 



cultural dynamics, etc 



yze tne dynamics 01 sucn processes m ( 

a y y y y y 0, y 0, y y 



l20j [. Even though 



in general the models studied in those works are simple ones, most of the results have been 
obtained via simulations. However, it is extremely relevant to have some form of analytical 
insight. 

. n this work we analyze a simple opinion formation model, analogous to the one studied 
in 2l| consisting of two parties, A and B, and an "intermediate" group /, that we call 
undecided agents. As in 2]], we assume that the supporters of parties A and B do not 
interact among them, but only through their interaction with the group /, convincing one 
of its members through a Sznajd-like rule similarly to what was discussed in [10], that is 
within a mean- field treatment. However, we don't consider that members of / can convince 
those of A or B, but instead we assume that there is a nonzero probability of a spontaneous 
change of opinion from / to the other two parties and viceversa: I ^ A, and I ^ B. We 
will see that this probability of spontaneous change of opinion (implying the existence of 



a social temperature 



s 



22,|23]) inhibits the possibility of reaching a consensus. Instead of 
consensus, we find that each party has some statistical density of supporters, and there is 
also a statistical stationary number of undecided (/) agents. 

Our aim is to write a master equation for this toy model, and study its behavior via 



van Kampen's f2-expansion approach 



24j . After determining if, in this case, the conditions 



for the validity of using such an approach are fulfilled, and exploiting it, we could obtain 
the macroscopic evolution equations for the density of supporters of A and B parties, as 
well as the Fokker-Planck equation governing the fluctuations around such deterministic 
or macroscopic behavior. The same approach also offers information about the typical 
relaxation behavior of small perturbations around the stationary macroscopic solutions. 
The outline of the paper is the following. In the next Section we present the model, and 
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apply van Kampen's f2 expansion approach in order to obtain the macroscopic equation and 
the Fokker-Planck equation governing the fluctuations around the macroscopic behavior. 
In Section 3 we analyze the behavior of the fluctuations through the study of their mean 
values and correlations, and discuss the relaxation time of small perturbations. In Section 
4 we present some typical results and finally, in Section 5, some general conclusions are 
summarized. 

II. THE MODEL AND THE APPROACH 
A. Description of the model 

We consider a system composed of three different groups of agents 

> supporters of the A party, indicated by Na, 

> supporters of the B party, indicated by Nb, 

> undecided ones, indicated by TV/. 

The interactions we are going to consider are only between A and I, and B and I. That 
means that we do not include direct interactions among A and B. The different contributions 
that we include are 

• spontaneous transitions A — > 7, occurring with a rate ol\ Na] 

• spontaneous transitions I — > A, occurring with a rate a 2 Ni; 

• spontaneous transitions B — > /, occurring with a rate a 3 N B ; 

• spontaneous transitions I — > B, occurring with a rate a^Nj; 

• convincing rule A + I — > 2 A, occurring with rate ^N^Nj; 

• convincing rule B + I — > 2 B, occurring with rate ^-NbNj. 

As indicated above, here iVj is the number of agents supporting the party or group "i" (with 
% = A,B,I). We have the constraint iV^ + iV^ + Nj — N, where N is the total number of 
agents. Such a constraint implies that, for fixed N, there are only two independent variables 
N A and N B . By using this constraint, the rates indicated above associated to processes 
involving Nj, could be written replacing Nj — (N — N A — N B ). 

With the above indicated interactions and rates, the master equation for the probability 
P(Na, N b , t) of having populations Na and Nb at time t (due we have had populations N^ 
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and N B at an initial time t a (< t)), may be written as 

^- t P{N A , N B , t) = at(N A + 1)P(N A + 1, N B , t) + a 3 (N B + 1)P(N A , N B + 1, t) 
+ a 2 (N-N A -N B + l)P(N A -l,N B ,t) 
+ a 4 (N -N A -N B + 1)P(N A , N B - 1, t) 
+ ^(N A -l)(N-N A -N B + l)P(N A -l,N B ,t) 

+ ^(N B -l)(N-N A -N B + l)P(N A ,N B -l,t) 
- a x N A + a 3 N B + a 2 (N - N A - N B ) 

+a 4 (N — N A — N B + 1)] P(N A , N B , t). (1) 

This is the model master equation to which we will apply van Kampen's approach 



B. Van Kampen's expansion 

In order to apply van Kampen's approach, as discussed in j^, we identify the large 
parameter Q with N (assuming N ^> 1); and define the following separation of the N^s into 
a macroscopic part of size f2, and a fluctuational part of size ^2, 

n a = m A (t) + nk A (t), 

n b = n* B (t) + n*£ B (t), (2) 

and define the density p — |£ (in our case p = 1). We also define the "step operators" 

Elf(N { ) = f(N l + l), 

ErVTO = f(N { - 1), 

with /(iVj) an arbitrary function. Using the forms indicated in Eqs. (J2J), in the limit of 
Q ^> 1, the step operators adopt the differential form 24] 

+1 ( 1 \ * d 1 / 1 \ d 2 

with % = A, B. Transforming from the old variables (N A ,N B ) to the new ones (£a,£b), we 
have the relations 

P(N A ,N B ,t) - U(U^B,t), (4) 
&—P(N A ,N B ,t) = —U(U,CbJ). (5) 



Putting everything together, and considering contributions up to order f^, yields the 
following two coupled differential equations for the macroscopic behavior 
d 



B 



(6) 
(7) 



It can be proved that the last set of equations has a unique (physically sound) stationary 
solution, i.e. a unique attractor 



V A (t -> oo) 

This is the main condition to validate the application of van Kampen's f2-expansion approach 



The following order, that is Q°, yields the Fokker-Planck equation (FPE) governing the 
fluctuations around the macroscopic behavior. It is given by 



d_ 

dt 



+ 



_d_ 

d 



1 

+ 2 



(a^ A + (« 2 + Pi^a)(U + 60 - PiU(p -*a- n(&, £s, *) 
(asCs + («4 + &*b)(£a + £b) - Mb(p -*a- *b)) n(a, 6?, *) 

«9 2 



«i^a + (a 2 + &*a)(p - * a - *b) 



1 

2 L 



d 2 



(8) 



As is well known for this approach |2J] , the solution of this FPE will have a Gaussian form 
determined by the first and second moments of the fluctuations. Hence, in the next section 
we analyze the equations governing those quantities. 



III. BEHAVIOR OF FLUCTUATIONS 



From the FPE indicated above (Eq. (jHJ)), it is possible to obtain equations for the 
mean value of the fluctuations as well as for the correlations of those fluctuations. For the 
fluctuations, (£a(£)) = Va and (£b(£)) = T]b, we have 



d . , 



d_ 

It 



VB(t) 



an + a 2 + Pi(2^ A + Vb) - Pip\ Va - (a 2 + Pi^a)Vb (9) 
a 3 + a 4 + f3 2 (^ A + 2^ B ) - f3 2 p r] B - (a 4 + (3 2 ^ b )va- (10) 
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Calling a a = (£, A (t) 2 ), &b = (£b(^) 2 ), an d &ab = (6a(*)£b(£)), we obtain for the correla- 
tion of fluctuations 

jfA{t) = -2ol x <t a - 2[a 2 + /3i9 A ] [a A + a AB \ + 2(3 1 a A \p - * A - * B ] 

+[ai* A + (o!2 + /9i*A)(p-*A-*B)], (11) 

Jfsif) = -2a 3 a B - 2[a 4 + Aj^b] [o^b + cr B ] + 2/? 2 <xb[p - ^ A - V B } 

+[a 3 V B + (a 4 + p 2 y B )(p-y A -y B )}, (12) 
-^ AB (t) = -[ati + a 3 ]o-Ai? - [«2 + A^a][ctab + 0b] 

-[04 + /^b]^ + v A b] + [p - ^ A - VbWPx + H°ab- (13) 



A. Reference state: symmetric case 

Here we particularize the above indicated equations to the symmetrical case, i.e. the case 
when = \E^. Hence, we adopt 

ot\ = «3 = a, «2 = «4 = a', 

and 

In such a case, the macroscopic equations © and (J7J) take the form 
d 



[a + a' - Z?]^ - - - a'^B + a' 



fit 

j t ^ B {t) = -[a + a' - p\9 B - - (3^ A ^ B - a'V A + a'. 



(14) 
(15) 



In order to make more explicit the solution of these equations, we work with the auxiliary 
variables S = ^ A + ^/ B and A = ty A — ty B , and use p = 1. Hence, the last equations 
transform now into 



a + 2a'-p £-/3£ 2 + 2o / 



d_ 

dt 



A(t) 



a — 13 



(16) 
(17) 



In the long time limit, t — > oo, we found on one hand 



A si = 0, 



implying ty A = \I/g, while on the other hand 



= pY? + 



a + 2a' - (3 



E - 2a'. 



This polynomial has two roots, but only one is physically sound, namely 



a + 2a' - (3 



-1 + 




wp 



[a + 2a' - p] 2 



(18) 



2^ • 



yielding = = vp** 

In a similar way, we can also simplify the equations for r] A and r) B , calling S(t) = r]A + T] B 
and D(t) = r] A — Vb- The corresponding equations are then rewritten as 



a + 2a' + 2/3(^ A + *b) - P 



d 
dt 



D(t) 



a + P{V a + Vb) ~ P 



D-P 



S. 



s. 



(19) 
(20) 



while for the correlation of the fluctuations we have 
d 



dt 

d 
It 



<TA(t) 
<T B (t) 



d . . 
Jt a AB {t) 



2aa A - 2[a' + pty A ] [a A + a AB ] + 2/3 [1 - V A - ^ B ]a A 
+ [a* A + (<*' + P^ A ){1 - V A ~ *b)] , 
2aa B -2[a' + P® B ] [a AB + a B ] + 2p[l -V A - V B ]a B 
+ [a^ B + (a' + P$ B )(1 - *a - *b)] , 
2aa j4i j - [a' + p^ A ] W AB + <?b] 
-[a' + pV B }[a A B + cr A ] + 2/?[l -V A - ^b]ctab- 



(21) 
(22) 
(23) 



Equations (fTH|) and (f2T)|) show that, in the asymptotic limit, i.e. for t — > oo, both, S = 
and -D = 0, implying that r/^* = ryg = 0. However, also in the general (non symmetric) case 
we expect to find r/ A = r/g = 0. In addition, from Eqs. (J2"T|) . (|2*2*j) and (J2HJ), it is clear that 
in general we obtain, again for t — > oo, that af ^ (i = A, B, AB). 

As we have seen, in the symmetric case we have ^ A = \l/g = hence it is clear that 
a A (t) and a B (t) behave in a similar way. And in particular a a = cr B = erf. In order to 
analyze the typical time for return to the stationary situation under small perturbations, we 
assume small perturbations of the form of « erf +5<7j(i) (i = A, B) and a AB ~ a AB +5<Ji(t), 



and fix ty s A 



\]>g = We find again that both Sa A (t) and 8a B (t) behave in the same way, 



and this help us to reduce the number of equations for the decay of correlations. Hence, we 
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can put Sa^it) = 5cr s (t) = 8a (t). The system driving the correlations becomes 
d 



dt Sa (t) 
j t 5a AB (t) 



-2 
-2 



a+a -p+3 pv* 
a + a' - (3 + 3 pV* 



5a Q - 
So ab - 2 



a' + /3V? 
a' + pm s l 



5a 



AB 



Sa Q . 



Clearly, 5af = 5a^ b = 0. After some algebraic steps we obtain 

5a (t) ~ 5or o (0) exp -2 [a + 2/3*f - /3] t 
5o-^ B (t) - *o-ab(0) exp [-2 [a + 2p% t - p] t 



(24) 
(25) 

(26) 
(27) 



These results indicate that, for the symmetrical case, the typical relaxation time is given by 



r relax = l -[a + 2pm s t -P}-\ 



(28) 



B. Beyond the symmetric case 

Let us call a a , a' Q and P a to the parameter's values corresponding to the symmetric case. 
We consider now the following cases where we vary the parameters 

Pi = Po, fo = Po + A/3, 
cti = a Q , «3 = a Q + Aa, 
«2 = Oi' , «4 = a' D + Aa'. 

We will vary only one of these parameters, while keeping the rest fixed. In the following 
section we present the results (mainly numerical) corresponding to those different cases. 



IV. RESULTS 

As indicated above, the macroscopic equations (Eqs. © and (|7|)) have a unique attractor, 
indicating that it is adequate to apply van Kampen's expansion approach. In this section 
we will present some results corresponding to symmetric and asymmetric situations, that 
show some typical behavior to be expected from the model and the approximation method. 
In what follows, al parameters are measured in arbitrary units. 

In Fig. we show the evolution of ^A{t) and ^/_B(t), the macroscopic solutions, indicating 
some trajectories towards the attractor: (a) for a symmetric, and (b) an asymmetric case. 
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FIG. 1: Evolution of the macroscopic solutions (Eqs. ([fi|7)l ). Case (a) corresponds to trajectories 
towards a symmetric solution (i.e. with = with parameters a\ = = 1, a.% = = 3, 

and /3i = p% = 2. Case (b) corresponds to trajectories towards an asymmetric solution (i.e. with 
7^ tyg), with parameters oti = 1, «3 = 5, 02 = «4 = 3, and Pi = P2 = 2. 

It is worth recalling that ^a and ^b are the density of supporters of party A and party 
B, respectively. During the evolution towards the attractor, starting from arbitrary initial 
conditions, we observe the possibility of a marked initial increase of the macroscopic density 
for one of the parties, follow by a marked reduction, or other situations showing only a 
decrease of an initial high density. Such cases indicate the need of taking with care the 
results of surveys and polls during, say, an electoral process. It is possible that an impressive 
initial increase in the support of a party can be followed for an also impressive decay of such 
a support. 

We remark that, due to the symmetry of the problem, it is equivalent to vary the set of 
parameters (a 3 , a 4 , /3 2 ) or the set («i, a 2 , Pi). Also worth remarking is That in both panels 
of Fig ^ the sum of ^a and is always < 1, so verifying that there is always a finite 
fraction of undecided agents. 

In Fig. 121 we depict the dependence of the stationary macroscopic solutions on different 
parameters of the system. On Fig. Efa) the dependence on is represented. It is apparent 
that for «3 < ati, we have ^/^ < \E'^, while for a% > ai, we find the inverse situation. 
Clearly, \1>^ = ^ when 03 = 1(= «i), as it corresponds to the symmetric case. Similarly, 
in Figs. El^b) andEfc) we see the dependence of the stationary macroscopic solutions on the 
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FIG. 2: Dependence of the stationary macroscopic solutions on different system parameters: (a) 
on a3, the rest of parameters are ot\ = oli = ay = 1, and (3\ = /?2 = 1- (b) on 04, the rest 
of parameters are a\ = = «3 = 1, and (5\ = = 1- ( c ) on /^2 5 the rest of parameters are 
ay = «2 = CK3 = ay = 1, and /?i = 1. In all three cases, the continuous line corresponds to 
while is indicated by the dotted line. 

parameters ay and 02, respectively. Also in these cases we observe similar behavior as in 
the previous one, when varying the indicated parameters. The parameters 03 or 04 (and 
similarly for a± or a 2 ) correspond to s pon taneous changes of opinion, and may be related 
to a kind of social temperature However, also (3\ and $2 are affected by such 

a temperature. So, the variation of these parameters in Fig. |21 correspond to changes in 
the social temperature, changes that could be attributed, in a period of time preceding an 
election, to increase in the level of discussions as well as the amount of propaganda. 

In Fig. El we depict the dependence of the stationary correlation functions for the fluctu- 
ations o~i (with % = 1, 2, corresponding to the projection of <ja,b,ab on the principal axes), on 
different systems' parameters. In Fig. Efa) the dependence on is represented, and simi- 
larly in Figs. Efb) and^c), the dependence on the parameters a 4 and /5 2 , respectively. We 
observe that, as the parameters are varied (that, in the case of 0-3 and ay, and as indicated 
above, could be associated to a variation of the social temperature) a tendency inversion 
could arise. This indicates that the dispersion of the probability distribution could change 
with a variation of the social temperature. This is again a warning for taking with some 
care the results of surveys and polls previous to an electoral process. 
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FIG. 3: Dependence of the stationary correlation functions <7j (with i = 1,2) corresponding to the 
projection of cta,b,ab on the principal axes, on different parameters of the system: (a) on 03, the 
other parameters are Qi = 02 = «3 = 1, and (3\ = /?2 = 1- (b) on a^, the other parameters are 
«i = OL2 = «3 = 1, and f3\ = fa = 1. (b) on the other parameters are a\ = Q2 = «3 = 04 = 1, 
and /3i = 1. 

Figure |H shows the stationary (Gaussian) probability distribution (pdf) II^, £s) s * pro- 
jected on the original (iV^, iV^) plane. We show three cases: on the left a symmetrical case, 
the central one corresponds to an asymmetrical situation with a population of N = 100, 
and on the right the same asymmetrical situation but with a population of iV = 1000. This 
last case clearly shows the influence of the population number in reducing the dispersion 
(as the population increases). We can use this pdf in order to estimate the probability pi 
(i = A,B), of winning for one or the other party. It corresponds to the volume of the dis- 
tribution remaining above, or below, the bisectrix N^/N = Nb/N. In the symmetrical case, 
as is obvious, we obtain Pa = Pb = 0.5 (or 50%), while in the asymmetrical case we found 
p B = 0.257 (or 25.7%) and p B = 0.015 (or 1.5%) for N = 100 and N = 1000, respectively. 
These results indicate that, for an asymmetrical situation like the one indicated here, we 
have a non zero probability that the minority party could, due to a fluctuation during the 
voting day, win the election. However, in agreement with intuition, as far as N ^> 1, and 
the stationary macroscopic solution departs from the symmetric case, such a probability Pi 
reduces proportionally to iV -1 jg^j. 

In Fig. on the left, we show a typical result for the time evolution of the macroscopic 
solution towards an asymmetric stationary case. In the same figure, in the central part 
we find the associated time evolution of the correlation functions for the fluctuations, crj 



0.6 
0.5 
0.4 

; tT 0.3 
0.2 
0.1 
0.0 
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FIG. 4: Stationary, Gaussian, probability distribution II (£4, £,b) s1 projected on the original 
(Na,Nb) plane. On the let we have a symmetrical case with ot\ = 03 = 2, «2 = «4 = 1> 
Pi = P2 = 2, and the population is N = 100. The central plot shows an asymmetrical case, with 
a\ = 2 and 03 = 2.5, while ai = 04 = 1, /3i = /32 = 2, and the population is N = 100. On the 
right we have the same asymmetrical case as before, but now N = 1000, showing the dispersion's 
reduction of the Gaussian distribution. 

(with i = 1,2) corresponding to the projection of (Ja,b,ab 011 the principal axes, while on 
the right we show the evolution of the angle between the principal axes and the figure axes. 
The temporal reentrance effect that has been observed in other studies exploiting the van 



Kampen's approach 



24. 



is apparent. This is a new warning, indicating the need to take 
with some care the results of surveys and polls during an electoral process. 

In Fig. IHlwe depict the dependence of the dominant (or relevant) relaxation time, that 
is the slowest of the three relaxation times, on different parameters of the system. On 
the left, we show a symmetrical case where the different lines represent the dependence 
respect to variation of: a.\ = a 3 indicated by a continuous line; a 2 = a 4 indicated by dotted 
line; (3\ = $2 indicated by dashed line. The strong dependence of the relaxation time on 
a = a 1 = «3 is apparent (in order to be represented in the same scale, the other two cases 
are multiplied by 3 or 10, respectively). This means that changes in the social temperature 
that, as discussed before, induce changes in a(— ct\ = a^), could significatively change the 
dominant relaxation time. On the right we show an asymmetrical case where, as before, the 
different lines represent the dependence respect to variation of: a\, indicated by a continuous 
line; a^, indicated by a dotted line; and /3i, indicated by dashed line. It is worth remarking 
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FIG. 5: On the left, we have the time evolution of the macroscopic solutions ^A(t) and The 
parameter values are ol\ = 1, = 5, a% = «4 = 3, Pi = 02 = 2. The stars indicate the position 
where the maxima that appear in the next panel occurs. Central part, time evolution of the 
correlation functions <7j (with i = 1,2) corresponding to the projection of (JAB AB on the principal 
axes. On the right, the angle between the principal axes and the figure axes. The parameters are 
a.\ = 1, «3 = 5, a?2 = «4 = 3, and 0\ = 0% = 2. 

that, when all the the parameters (a%, a 2 and f3{) are equal to 1, we see that the relaxation 
time is the same. On the left figure, this is shown in the inset. In the asymmetrical case, 
the behavior is of the same order for the variation of the three parameters. However, the 
comment about the effect of changes in the social temperature remain valid. 



V. CONCLUSIONS 



We have studied a simple opinion formation model (that is a toy model), analogous to 



the one studied in 



2l| . It consists of two parties, A and B, and an intermediate group /, 



that we call undecided agents. It was assumed that the supporters of parties A and B do 
not interact among them, but only through their interaction with the group /, convincing 
its members through a mean- field treatment; that members of / are not able to convince 
those of A or B, but instead we consider a nonzero probability of a spontaneous change of 
opinion from I to the other two parties and viceversa. It is this possibility of spontaneous 
change of opinion that inhibits the possibility of reaching a consensus, and yields that each 
party has some statistical density of supporters, as well as a statistical stationary number 
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0-1 1 ' 1 ■ 1 ' 1 ■ 1 1 0,4-1 1 , 1 , 1 , 1 , 1 1 

0.0 0,5 1,0 1,5 2,0 o,0 0,5 1,0 1,5 2,0 

FIG. 6: Dependence of the dominant relaxation time on different system parameters. On the left, 
symmetrical case: continuous line varying a\ = 03, dotted line varying «2 = «4> and dashed line 
varying f3\ = /?2- In order to compare all three, the dotted line was multiplied by 3, while the 
dashed one by 10. The inset shows, now on the same scale, the crossing of the lines at the point 
where all the parameters are equal to 1. On the right, asymmetrical case: continuous line varying 
«i, dotted line varying 0*2, and dashed line varying (3\. In all cases, the parameters that remain 
constant are all = 1. 

of undecided agents. 

Starting from the master equation for this toy model, the van Kampen's f2-expansion 
approach |24| was exploited in order to obtain the macroscopic evolution equations for the 
density of supporters of A and B parties, as well as the Fokker-Planck equation governing 
the fluctuations around such a macroscopic behavior. Through this same approach infor- 
mation about the typical relaxation behavior of small perturbations around the stationary 
macroscopic solutions was obtained. 

The results indicate that one needs to take with care the results of social surveys and 
polls in the months preceding an electoral process. As we have found, it is possible that an 
impressive initial increase in the support of a party can be followed for an also impressive 
decay of such a support. The dependence of the macroscopic solutions as well as the corre- 
lation of the fluctuations on the model parameters, variation in a^, or j3 2 (that, due to 
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the symmetry of the model are similar to varying ai, a 2 or fc) was also analyzed. As the 
parameters correspond to spontaneous change of opinion, or $ to convincing capacity, 
and it is possible to assume that have an "activation-like structure" , we can argue that this 
could be related to changes in the social temperature, and that such a temperature could 
be varied, for instance, in a period near elections when the level of discussion as well as the 
amount of propaganda increases. 

We have also analyzed the probability that, due to a fluctuation, the minority party 
could win a loose election, and that such a probability behaves inversely to iV (the popu- 
lation number). Also analyzing the temporal behavior of the fluctuations some "tendency 
inversion" indicating that, an initial increase of the dispersion could be reduced as time 
elapses was found. 

We have also analyzed the relaxation of small perturbations near the stationary state, 
and the dependence of the typical relaxation times on the system parameters was obtained. 
This could shead some light on the social response to small perturbations like an increase of 
propaganda, or dissemination of information about some "negative" aspects of a candidate, 
etc. However, such an analysis is only valid near the macroscopic stationary state, but looses 
its validity for a very large perturbation. For instance, a situation like the one lived in Spain 
during the last elections (the terrorist attack in Madrid on March 11, 2003, just four days 
before the election day), clearly was a very large perturbation that cannot be described by 
this simplified approach. 

Finally, it is worth to comment on the effect of including a direct interaction between 
both parties A and B. As long as the direct interaction parameter remains small, the monos- 
tability will persist, and the analysis, with small variations will remain valid. However, as 
the interaction parameter overcomes some threshold value, a transition towards a bistability 
situation arise, invalidating the exploitation of the van Kampen's f2-expansion approach. 
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